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Abstract. Shock acceleration is an ubiquitous phenomenon in astrophysical plasmas. Plasma waves 
and their associated instabilities (e.g., Buneman, Weibel and other two-stream instabilities) created 
in collisionless shocks are responsible for particle (electron, positron, and ion) acceleration. Using 
a 3-D relativistic electromagnetic particle (REMP) code, we have investigated particle acceleration 
associated with a relativistic jet front propagating into an ambient plasma with and without initial 
magnetic fields. We find small differences in the results for no ambient and modest ambient 
magnetic fields. Simulations show that the Weibel instability created in the collisionless shock front 
accelerates jet and ambient particles both perpendicular and parallel to the jet propagation direction. 
The non-linear fluctuation amplitudes of densities, currents, electric, and magnetic fields in the 
electron-positron shock are larger than those found in the electron-ion shock at the same simulation 
time. This comes from the fact that both electrons and positrons contribute to generation of the 
Weibel instability. While some Fermi acceleration may occur at the jet front, the majority of electron 
and positron acceleration takes place behind the jet front and cannot be characterized as Fermi 
acceleration. The simulation results show that the Weibel instability is responsible for generating 
and amplifying nonuniform, small-scale (mainly transverse) magnetic fields which contribute to 
the electron's (positron's) transverse deflection behind the jet head. This small scale magnetic field 
structure is appropriate to the generation of "jitter" radiation from deflected electrons (positrons) 
as opposed to synchrotron radiation. The jitter radiation has different properties than synchrotron 
radiation calculated assuming a uniform magnetic field. The jitter radiation resulting from small 
scale magnetic field structures may be important for understanding the complex time structure 
and spectral evolution observed in gamma-ray bursts and other astrophysical sources containing 
relativistic jets and relativistic collisionless shocks. 



INTRODUCTION 



Nonthermal radiation observed from astrophysical systems containing relativistic jets 
and shocks, e.g., active galactic nuclei (AGNs), gamma-ray bursts (GRBs), and Galactic 
microquasar systems usually has power-law emission spectra. In most of these systems, 



the emission is thought to be generated by accelerated electrons through the synchrotron 
and/or inverse Compton mechanisms. Radiation from these systems is observed from 
the radio through the gamma-ray region. Radiation in optical and higher frequencies 
typically requires particle acceleration in order to counter radiative losses. It has been 
proposed that the needed particle acceleration occurs in shocks produced by differences 
in flow speed. 

Particle-in-cell (PIC) simulations can shed light on the physical mechanism of particle 
acceleration that occurs in the complicated dynamics within relativistic shocks. Recent 
PIC simulations using injected relativistic electron-ion jets show that acceleration oc- 
curs within the downstream jet, rather than by the scattering of particles back and forth 
across the shock as in Fermi acceleration D2L uL la> 120 > and Silva et al. 111(11 have presented 
simulations of the collision of two inter-penetrating electron-positron plasma shells as a 
model of an astrophysical collisionless shock. In the electron-positron simulations per- 
formed with counter-streaming jets 111 Of . shock dynamics involving the propagating jet 
head (where Fermi acceleration may take place) was not investigated. In general, these 
independent simulations have confirmed that relativistic jets excite the Weibel instability 
I llll . The Weibel instability generates current filaments and associated magnetic fields 
1 2D, and accelerates electrons 1 1 QL I2L I3L I8L 1411 . 

In this paper we present new simulation results of particle acceleration and magnetic 
field generation for relativistic electron-positron and electron-ion shocks using 3-D rel- 
ativistic electromagnetic particle-in-cell (REMP) simulations. In our new simulations, 
electron-positron and electron-ion relativistic jets with Lorentz factor, 7=5 (corre- 
sponds to 2.5 MeV) is injected into electron-positron and electron-ion plasmas in order 
to study the dynamics of a relativistic collisionless shock both with and without an initial 
ambient magnetic field. 



SIMULATION SETUP AND RESULTS 

Four simulations were performed using an 85 x 85 x 320 grid with a total of 180 million 
particles (27 particles/cell/species for the ambient plasma) and an electron skin depth, 
^ce = cj (Ope = 9.6A, where C0p e = (4-7te 2 n e /m e ) 1//2 is the electron plasma frequency and 
A is the grid size. In all simulations jets are injected at z = 25 A in the positive z direction. 
In all simulations radiating boundary conditions were used on the planes at z = 0,z max . 
Periodic boundary conditions were used on all other boundaries J1D. 

In two simulations an electron-positron jet is injected into a magnetized and un- 
magnetized electron-positron ambient plasma and in two simulations an electron-ion 
jet is injected into a magnetized and unmagnetized electron-ion ambient plasma. The 
choice of parameters and simulations allows comparison with previous simulations 
[10,2|iSIi. 

The electron number density of the jet is 0.741/ib, where is the density of ambient 
(background) electrons. The average jet velocity is v; = 0.9798c, and the corresponding 
Lorentz factor is 5. The jet is cold (vj th = vj^ = 0.01c and vj th = 0.0022c in the 
laboratory frame). The ambient and jet electron-positron plasma has mass ratio m p / m e = 
m e +/m e - = 1 (mi/m e = 20). The electron and ion thermal velocities in the ambient 



plasma are v^ h = 0.1c and v[ h = 0.022, respectively, where c is the speed of light. 
The time step t = 0.013/c0p e , the ratio (Op e /Q. e = H-5, and the Alfven speed (for 
electrons) VAe = (&e/(Dpe)c = 8.66 x 10~ 2 c. With the speed of an Alfven wave given 
by v A = Wl/il +Vl/c 2 )] 1 / 2 = 6.10 x 10~ 2 c where V A = [B 2 /4n(n s m s + n p m p )] 1 / 2 = 
6. 12 x 10~ 2 c, the Alfven Mach number Ma = vj/va = 16.0. With a magnetosonic speed 
v ms = (v 2 h + vjj 1 / 2 = 0.132c the magnetosonic Mach number M ms = Vj/v ms = 7.406. 
At least approximately the appropriate relativistic Mach numbers multiply these values 
by the Lorentz factor. Thus, in an MHD approximation we are dealing with a high Mach 
number shock with yM >> 1. The gyroradius of ambient electrons and positrons with 
v_l = vth = 0. lc is 1 1 . 1 A = 1 . 154A ce (for ambient ions: 49. 6A = 5.16A ce ). All the Mach 
numbers with electron-ion jets are approximately increased by ^Jm\Jrn e = \/2ti = A.A1. 




FIGURE 1. 2D images in the x — y plane at z — 230A for a fiat jet injected into a unmagnetized ((a) 
and (b)) and magnetized ((c) and (d)) ambient medium shown at t = 28.8/©p e . Colors indicate the z- 
component of the current density (J z ) ((a) and (c): electron-positron, (b) and (d): electron-ion) (peaks: (a) 
±54.7, (b) -123.6, (c) ±58.0, (d) -135.5 with B x ,B y indicated by the arrows. 
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FIGURE 2. Velocity distributions of jet electrons four cases as Fig. 1 at t — 28.8/a)pe ((a): electro- 
positron, unmagnetized, (b) electron-ion, unmagnetized, (c)electro-positron, magnetized, and (d) electron- 
ion, magnetized). The blue and red curves show the distributions of injected and shocked jet electrons. Jet 
electrons are binned as a function of yvii , where y= (1 — (vjj + v'jj/c 2 ) -1 / 2 . 



Current filaments resulting from development of the Weibel instability behind the jet 
front (at z = 230A) are shown in Figure 1 at time t = 28.8/(Op e for four different cases. 
The upper two panels (a and b) show for unmagnetized ambient plasmas and the lower 
two panels (c and d) are for magnetized ambient plasmas. In cases (a and c) an electron- 
positron jet is injected in an electron-positron ambient plasma. An electron-ion jet is 
injected in an electron-ion ambient plasma (b an d). The maximum values of J z are (a) 
±54.7, (b)— 123.6, (c) ±58.0, and (d) —135.5. The electro-positron jets show electron 
(negative) and positron (positive) current filamentations, since both species contribute 
to the Weibel instability. On the other hind, at this simulation time mainly electron jets 
generate (negative) current filamentations as shown in Figs. 1 b and Id. The effect of 
weak ambient magnetic fields increase the maximum values of current filamentations 



slightly. 

The electrons are deflected by the transverse magnetic fields (B x ,B y ) via the Lorentz 
force: — e(v x B), generated by current filaments (7 Z ), which in turn enhance the trans- 
verse magnetic fields lfi"lL l7ll. The complicated filamented structures resulting from the 
Weibel instability have diameters on the order of the electron skin depth (A ce = 9.6A). 

This is in good agreement with the prediction of X ~ 2 1//4 cy t j 1 //2 / C0p e ~ 1 . 1 88 A ce = 11 .4A 
0]. Here, y t h ~ 1 is a thermal Lorentz factor. However, in the electron-positron jets the 
current filaments are coalesced in the transverse direction, and this shows the nonlin- 
ear evolution. The longitudinal current (7 Z ) in the electron-positron jets (a and c) shows 
significantly more transverse variation than in the electron-ion jets (b and d). 

The acceleration of electrons has been reported in previous work iTToL 0, S S S 01- 
Figure 2 shows that the cold jet electrons are accelerated and decelerated. As expected, at 
this time jet electrons in the electron-positron jets are thermalized more strongly than in 
the electron-ion jets. The blue curves in Figs. 2b and 2d is close to the initial distribution 
of injected jet electrons. We also see that the kinetic energy (parallel velocity vn ~ vj) of 
the jet electrons is transferred to the perpendicular velocity via the electric and magnetic 
fields generated by the Weibel instability [9]. The strongest transverse acceleration of jet 
electrons accompanies the strongest deceleration of electron flow and occurs between 
z/A = 210 — 240. The transverse acceleration in the electron-positron jets is over four 
times that in the electron-ion simulations. The strongest acceleration takes place around 
the maximum amplitude of perturbations due to the Weibel instability at z/A ~ 220 as 
seen in Figs la and lc. 

SUMMARY AND DISCUSSION 

We have performed self-consistent, three-dimensional relativistic particle simulations 
of relativistic electron-positron and electron-ion jets propagating into magnetized and 
unmagnetized electron-positron and electron-ion ambient plasmas. The main acceler- 
ation of electrons takes place in the downstream region. Processes in the relativistic 
collisionless shock are dominated by structures produced by the Weibel instability. This 
instability is excited in the downstream region behind the jet head, where electron den- 
sity perturbations lead to the formation of current filaments [9] . The nonuniform electric 
field and magnetic field structures associated with these current filaments thermalize 
the jet electrons and positrons, while accelerating the ambient electrons and positrons, 
and accelerating (heating) the jet and ambient electrons and positrons in the transverse 
direction. 

Other simulations with different skin depths and plasma frequencies show that the 
growth and structure of current filaments generated by the Weibel instability scale with 
the plasma frequency and the skin depth |9j]. An additional simulation in which an 
electron-ion jet is injected into a ambient plasma with perpendicular magnetic field 
shows magnetic reconnection due to the generation of an antiparallel magnetic field 
generated by bending of jet electron trajectories Js|]. 

This small scale magnetic field structure generated by the Weibel instability is ap- 
propriate to the generation of "jitter" radiation from deflected electrons (positrons) Hal. 



The jitter radiation resulting from small scale magnetic field structures needs to be cal- 
culated from the trajectories of electrons and positrons (ions) which requires extensive 
computational resources. This investigation will provide an idea for understanding the 
complex time structure and spectral evolution observed in gamma-ray bursts and other 
astrophysical sources containing relativistic jets and relativistic collisionless shocks. 

The fundamental characteristics of relativistic shocks are essential for a proper un- 
derstanding of the prompt gamma-ray and afterglow emission in gamma-ray bursts, and 
also to an understanding of the particle reacceleration processes and emission from the 
shocked regions in relativistic AGN jets. Since the shock dynamics is complex and sub- 
tle, more comprehensive studies using larger systems are required to better understand 
the acceleration of electrons, the generation of magnetic fields and the associated emis- 
sion. This further study will provide the insight into basic relativistic collisionless shock 
characteristics needed to provide a firm physical basis for modeling the emission from 
shocks in relativistic flows. 
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